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Abstract. The dipole modes of non-parabolic quantum dots are studied by means of their current and 
density patterns as well as with their local absorption distribution. The anticrossing of the so-called Bern- 
stein modes originates from the coupling with electron-hole excitations of the two Landau bands which 
are occupied at the corresponding magnetic fields. Non-quadratic terms in the potential cause an energy 
separation between bulk and edge current modes in the anticrossing region. On a local scale the fragmented 
peaks absorb energy in complementary spatial regions which evolve with the magnetic held. 

PACS. 73.21.B Electron states and collective excitations in multilayers, quantum wells, mesoscopic and 
nanoscale systems - 73.20.Mf Collective excitations (including excitons, polarons, plasmons and other 
charge-density excitations) 



1 Introduction 

The far infrared (FIR) spectroscopy has proved to be an 
invaluable tool for the physical characterization of semi- 
conductor quantum dots and other electronic nanostruc- 
tures [0,0]. Since its initial applications to 2D semiconduc- 
tor quantum dots one of the main motivations has been 
to identify a signal of the relative motion of the confined 
electrons. As is well known [3|, most commonly the con- 
finement seems to be parabolic, and in this limit the gen- 
eralized Kohn theorem assures that the only allowed ex- 
citations are the center-of-mass modes, at the frequencies 
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where B is the applied perpendicular magnetic field, Uq is 
the frequency associated with the external parabola and 
uj c = eB/c is the cyclotron frequency [Q. 

In Ref . |^| the fragmentation of the high-energy branch 
uj-\-{B) was measured in wires and dots and interpreted 
as an interaction with the cyclotron harmonics at ener- 
gies nuj c (n — 2,3). Since these interactions resemble the 
Bernstein modes of the electron gas ](| , the same name was 
used to label the corresponding excitations in these nanos- 
tructures. More recently, Krahne et al. (7) have measured 
the Bernstein modes in a 2D GaAs system with tunable 
electron density, varying from a continuous 2D gas to well 
separated dots. While the anticrossing exactly occurs at 
2uj c for the 2D gas, for finite quantum dots it lies below 
(between lu c and 2ui c ). 



It is our aim in this paper to provide a physical charac- 
terization of the below- 2u> c Bernstein modes of quantum 
dots by analysing the current and density distributions, as 
well as the local absorption patterns, associated with each 
particular peak of the dipole spectrum. For this purpose, 
we shall use the time-dependent local-spin-density approx- 
imation in a symmetry unrestricted formalism as devel- 
oped recently in Ref. || . A comparison with the electron- 
hole interband transitions will also shed light on the origin 
of the fragmentation of the uj + (B) branch. 



2 Theoretical approach 

As in Ref. || we assume the following confining potential 
1 



"confM = -LO 



to 
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with r the radial coordinate in the xy plane, where the 

— 1/2 

electronic motion occurs, £q = lo the confining length, 
u>o = 3.37 meV and a — 2.02 1CP 2 . The parameter a con- 
trols the degree of non-parabolicity although, as shown in 
Ref. H , the Bernstein fragmentation discussed below is ac- 
tually not much sensitive on its precise value. It is worth 
to mention that a similar potential including a quartic 
term was used by Ye and Zaremba to analyze the break- 
ing of Kohn's theorem in the context of a hydrodynamic 
approach ||. 

We assume the effective-mass Hamiltonian, applicable 
to GaAs nanostructures |l0|, and describe the electronic 
exchange and correlation effects within the Local-Spin- 
Density approximation (LSD A). This density-functional 
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approach has recently been used by several authors to de- 
scribe quantum dot properties. The reader is addressed 
to Refs. 0,|l2],[l3) for details on the method. In the con- 
text of the FIR absorption the time-dependent extension 
(TDLSDA) has also been shown to provide adequate re- 
sults 

The dynamical properties described below have been 
calculated by integrating the Kohn-Sham equations 



.dip iv (r,t) 

i — = h n [p,m\ tp in [T,t) 



(3) 



where 77 =|, J. labels the two spin components, while p = 
\ifil\ 2 indicates the total particle density 
|<^i|| 2 — IviJ 2 the total spin magnetiza- 



and m = V,,, 

tion. The self-consistent Hamiltonian in Eq. [3| reads 

1 / 6 \ 2 

h v [p, m] = - (^-iV + - A ( r )J + v COI a(r) 



dt 



, P(r') , SE X c[p,m] 
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Sr, I -z Vg pB-r 

dm I 
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with A(r) = B/2(—y,x) being the vector potential in 
the symmetric gauge, ExcIPt 11 ^] the exchange-correlation 
functional, and = ±1 for rj =] . j, respectively. The 
Zeeman term in Eq. ^ contains the gyromagnetic factor 
g* and the Bohr magneton ps = e/2m e c. 



3 Single peak analysis 

After calculating the ground state structure and in order 
to excite the collective dipole mode we perform a small 
rigid translation of the electronic cloud. This takes the 
system out of equilibrium and as a consequence it begins 
to oscillate. The expectation values of several observables 
O are subsequently recorded in time (0)(t) and frequency 
analysed to obtain the corresponding energy distributions. 
In the present analysis we have considered the dipole, local 
current and local density operators, i.e., 



D = ^2 Xi + yi , 

i 

J(r)=]T ~(Vi-V, 

i L 

p(r) =^<5(r-r i ) . 



-A(r,) 

c 



(5) 



Note that after the frequency analysis we obtain (0)(ui) 
which for the local signals provides a pattern of currents 
(j(r))(cj) and of local density (p(r)}(u) that can obviously 
be ascribed to the excitation at energy u>. At a fixed u>, the 
time evolution at every point r is simply given by a phase 
e -1 "*, or by sm(ut) and cos(wi) for the corresponding real 
transforms, thus permitting to monitor the variation of the 
current and density patterns for each excitation peak. In 



addition, the density patterns can be used to obtain the 
local absorption from | (p(r))(u)) |, as done in Ref. |r| for 
the characterization of the FIR absorption of triangular 
and square quantum dots. Theoretically, the local absorp- 
tion at R provides the energy absorbed by the system 
when a probe of the type £(r — R) is used, where £ is a 
highly peaked spatial modulation p7| . Note also that in 
Eq. we have used the gauge invariant current including 
explicitly the vector potential A(r). 



4 Results 

Figure 1 displays the FIR absorption for a dot with N = 
20 electrons in the confining potential (||) for different ver- 
tical magnetic fields. The Bernstein fragmentation of the 
high energy branch is conspicuous for 1.6T< B < 2.4 T. 
Actually, for several B's up to three peaks can be identi- 
fied in the higher branch, although in most instances this 
branch fragments into two dominant features. In this fig- 
ure the dashed line indicates the 2cu c values while the dot- 
ted lines show the particle-hole transitions with Al = +1, 
which by angular-momentum selection rules are the only 
ones that couple to the high-energy branch. All these tran- 
sitions are of interband character. The At = — 1 particle- 
hole transitions (of both intra- and interband character) 
only contribute to the low-energy branch and thus are 
not relevant for the present discussion. From this figure 
and the results that will be presented below, the frag- 
mentation of the FIR absorption, causing the Bernstein 
modes, can be understood as an effect of the coupling with 
particle-hole transitions, allowed by the non-parabolicity 
of the confining potential. Note that at B w 2 T there 
is a particle-hole excitation actually lying in between the 
fragmented Bernstein peaks, thus showing that the maxi- 
mal overlap of the interacting response with the interband 
transitions occurs below the 2lo c line. We stress that in 
spite of the overlap with oj-(B), the displayed transitions 
do not couple with the lower branch by angular momen- 
tum selection rules. 

The relevance of the different particle-hole transitions 
is more easily appreciated within the perturbative response 
formalism, as opposed to the real-time one (see, e.g., Ref. 
[p!s| ). Using the perturbative method, in Fig. 2 we com- 
pare the full absorption function at B = 1.8 T with those 
obtained by including only the hole states from the first 
and second Landau bands, respectively, shown in the left 
panel. The interacting responses in the restricted sub- 
spaces show collective peaks which approximately repro- 
duce the energies of the Bernstein excitations in the full 
response. Therefore, each of the Bernstein peaks can be 
thought of as arising from particle-hole excitations of dif- 
ferent bands, as it happens in the bulk limit for a mod- 
ulated 2D gas |iS|| . Note, however, that the picture of 
separate FIR excitations for each band is an approxi- 
mate one, since there are important interference effects 
between bands. With this interpretation we expect that 
when lo + (B) and 2uj c (B) intersect in a region where only 
one Landau band is occupied there will be no Bernstein 
fragmentation as in Fig. 1. Indeed we have checked this 
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by using a much lower value of the confinement, ujq = 
0.5 meV, which for B = 0.3 T has one (spin degenerate) 
occuppied Landau band and uj+(B) ~ 2u c (B) without 
any anticrossing. 

The preceding analysis is in qualitative agreement with 
the results for Raman modes in quantum wires of Steinebach 
et al. [^(| . These authors attribute the anticrossing to the 
large importance of transitions with An — 2, where n is 
the Landau level index, with respect to those with An = 1 . 
Realizing that An = 1 transitions are mainly from the 
highest occupied Landau band, while An — 2 also in- 
cludes the second highest Landau band, one arrives at a 
similar conclusion as the above one. 

We consider next the current distributions. The pure 
Kohn modes in parabolic dots correspond to rigid trans- 
lations of the electronic density and, therefore, to essen- 
tially uniform current distributions, the only variations 
being due to the density inhomogeneities. A similar result 
is found for the low energy branch in Fig. 1, indicating the 
Kohn-like character of this mode. 

The current patterns at a given time for the Bern- 
stein modes in the nonparabolic dot are shown in Fig. 3. 
The time evolution for each of these patterns is simply an 
anticlockwise rotation, in agreement with the Ai = +1 
character of the mode (the positive z axis is pointing to- 
wards the reader). The leftmost panel corresponds to the 
dominant peak at B = 1.4 T, which displays a current 
pattern basically uniform in the bulk of the dot, thus sim- 
ilar to the previously mentioned Kohn modes. The two 
intermediate panels correspond to the peaks at B = 2 T 
(at w « 0.52 and 0.54). They show an incipient separa- 
tion of bulk and edge current oscillations which is fur- 
ther developed at B = 2.4 T (two right panels). Note 
that the edge-current patterns contain a hole in the dot 
center. Therefore, the non-parabolicity of the potential in- 
duces a separation in energy of the higher-branch bulk and 
edge current modes, which are degenerate in the purely 
parabolic case. 

It is also interesting to look at the density variation 
patterns displayed in Fig. 4 for the B = 2.4 T case. As 
before, the time evolution of each pattern is a clockwise ro- 
tation for the lower-energy branch and anticlokwise for the 
upper modes. Red and blue indicate increment and decre- 
ment in the local density, respectively. Note that both the 
low-energy peak (fl) and the dominant peak of the higher 
branch (f2) exhibit a simple dipole pattern with two lobes. 
On the contrary, the highest mode (f3) displays an internal 
structure reflecting the existence of a node in the radial 
density and four regions with alternate phases. The inter- 
nal structure of f3 nicely correlates with the the current 
pattern of Fig. 3, since the regions of current convergence 
(divergence) correspond to an increase (decrease) of the 
local density. 

A time average of the oscillating density amplitudes 
highlights the regions of higher FIR absorption in the dot, 
as shown in Fig. 5. The different grey colours indicate the 
absorption strength as given by the local oscillation. We 
notice that at B = 1.4 T the Bernstein peaks absorb basi- 
cally in rings, with the lower peak having a more internal 



character. As the magnetic field is raised the lower peak 
tends to expand its absorption ring, while the higher one 
absorbs more in the inner region. As a result, at B w 2.6 T 
the two modes reverse their character, the upper mode be- 
coming more internal. Thus showing that the two Bern- 
stein peaks absorb energy in complemcmentary spatial re- 
gions. 



5 Conclusions 

We have analysed the fragmentation of the high-energy 
branch in non-parabolic dots by using a) a comparison 
with the allowed electron-hole transitions, b) the current 
and density variations associated with each peak, and c) 
the local absorption pattern of each mode. We conclude 
that the Bernstein fragmentation is a result of the coupling 
with electron- hole dipole transitions originating from each 
of the two occupied Landau bands. This coupling mani- 
fests in an energy separation of the bulk and edge current 
modes, with the latter one having a four-lobe structure in 
the oscillating density. The separation of the modes is also 
reflected in the local absorption, being more internal for 
the lower peak than for the higher one at small B's and 
reverting character when the magnetic field is increased. 

This work was supported by Grant No. PB98-0124 
from DGESeiC, Spain, the Research Fund of the Univer- 
sity of Iceland, and the Icelandic Natural Science Council. 
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Fig. 1. Evolution of the dipole strength for the 20-electron dot. The dashed line indicates ui = 2uj c and the dotted one the 
Ai = +1 electron-hole interband transitions. The labels identify the different peaks discussed in the text. 
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Fig. 2. Left panel displays the Kohn-Sham single particle energies as a function of angular momentum at a magnetic field 
B = 1.8 T. The dashed line separates occupied and non-occupied levels. Right panels show the response using the perturbative 
approach |lq). The dashed lines correspond to the restricted subspaces for hole states in Landau bands A and B, respectively. 



M. Valm- Rodriguez et al: Bernstein modes in quantum dots. . . 



6 - 



- 



-6 - 



///././ a ^ 

/ JJJJJJ,/ , 
/ JJJJJJ,/, 



/ /////// / 
/////// 



I 

6 



. , 1 / / //. 




, 1 1 


f3 




» 1 f 


. / / / / / /s . 




/ / / 


/ / // 


. / / / s s / / 




, ^ ^ / ' > 


, / / 


















_ ^/ /,,/// ^ . 








, ^/V s s / / / . 




*y S / * > 




. /////// , 




s // // 


/ / 


,////!« . 




, /// / 


1 , 


.//II, 




/ / / 


f - 











d2 




f2 


. / / / / / , , 




/ / / / / < 


/ / /////,. 




,///////, 


, / / //////, 




.////////, 


, / ///////, . 




,/////////, 


, ,////////, 




,/////////, 


, , ///////> , 




, / ////////, 


, /////// 1 , 




,////////. 


. , ////// / 




, * /////// 


, y / / / / / , 




✓ / / / / * 


r * S S 4> 







Fig. 3. Current density patterns for the high energy branch at different magnetic fields. Each label indicates the corresponding 
mode in Fig. f . A different absolute scale, varying in proportion to the dipole strength of the corresponding peak of Fig. f , has 
been used in each panel. 



Fig. 4. Pattern of local-density variations for the indicated modes of Fig. 1 which correspond to the three dipole peaks at 
B — 2.4 T. The gray color scale indicates the signal magnitude while the corresponding sign is superimposed in white. 



Fig. 5. Local absorption patterns for the Bernstein peaks at different magnetic fields. Black regions correspond to high absorp- 
tion, with an absolute scale proportional to the peak's dipole strength. 
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